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A novel algebra underlying integrable systems is shown to generate and unify a large class of 
quantum integrable models with given _R-matrix, through reductions of an ancestor Lax operator and 
its different realizations. Along with known discrete and field models a new class of inhomogeneous 
and impurity models are obtained. 
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on : The self-dual Yang-Mills equation with possible reductions has given a vivid unifying picture in classical integrable 
systems in 1 + 1 and + 1 dimensions [Q. However, in the quantum case not much has been achieved in this direction 
and there exists a genuine need for discovering some scheme, which would generate models of quantum integrable 
' systems (QIS) Q along with their Lax operators and i?-matrix and thus unify them. 

ON ' The significance of algebraic structures in describing physical consequences is well recognized. Like Lie algebras, 
their quantum deformations || also found to be of immense importance in physical models B- Q. In fact the idea 

£N| of quantum Lie algebra, which attracted enormous interest in recent years |7j- || has stemmed from the QIS and at 
the same time has made profound influence on the QIS itself flip]] - flfj}| . 

£NJ | Motivated by these facts and our experience jl2| , we find a novel Hopf algebra as a consequence of the integrability 
condition, which underlies integrable models with 2x2 Lax operators and the trigonometric i?-matrix. This is more 

y—i • general than the well known quantum Lie algebra and in contrast represents a deformed quadratic algebra (QdA), 



so called due to the appearance of generators in quadratic form in the defining algebraic relations. At the same time 
it unifies a large class of quantum integrable models by generating them in a systematic way through reductions of 
an ancestor model with explicit Lax operator realization. Note that the Lax operator together with the quantum 
O || i?-matrix define an integrable system completely, giving also all conserved quantities including the Hamiltonian of the 
model. 

The proposed algebra may be given by the simple relations 



X' [S ,3 .S' ± ] =±5*, \S+,S-} = (M+ sin(2aS 3 ) + M- cos(2aS 3 )) , [Af ± ,-1 = 0, (1) 

where are the central elements. We show that (Q) is not merely a modification of known U q (su(2)), but a QdA 
underlying an integrable ancestor model and in effect is dictated by the quantum Yang Baxter equation (QYBE) |Tij] 
RLL = LLR. The associated quantum i?(A)-matrix is the known trigonometric solution related to sine-Gordon (SG) 
Q , while the Lax operator may be taken as 

4 2sinaS+ + f^e*"** J ' 4 ' { ) 

with central to ([!]) relating = ±v / ±T(c^ l "c 2 ' ± cj~c^). The derivation of algebra (Q) follows from QYBE by 
inserting the explicit form (Q) and the i?-matrix and matching different powers of the spectral parameter £. 

Note that (0) is a Hopf algebra Jl8| and a generalization of U q (su{2)). However, unlike Lie algebras or their 
deformations, due to the presence of multiplicative operators M^, (^) becomes quantum-deformation of a QdA. 
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Since these operators have arbitrary eigenvalues including zeros, they can not be removed by scaling and therefore 
generically (|l]) is different from the known quantum algebra. Moreover different representations of M's generate new 
structure constants leading to a rich variety of deformed Lie algebras, which are related to different integrable systems. 
This fact becomes important for its present application. The appearance of QdA in basic integrable system should be 
rather expected, since the QYBE with i?-matrix having c-number elements is itself a QdA. The notion of QdA was 
introduced first by Sklyanin Jl9[ . 

The ancestor model can be constructed through representation of ([[]) in physical variables ( with [u,p] — i,) as |20[| 

S 3 = u, S+ = e- ip g(u), S~=g{ U y p . (3) 

This gives a novel exactly integrable quantum system generalizing lattice SG model and associated with the Lax 
operator (||). ft is evident that for hermitian g(u) only one gets S~ = (S + Y . We show below that through various 
realizations of the single object g(u) in (||) the ancestor model generates a whole class of integrable models. Their 
Lax operators are derived from (^|), while the i?-matrices are simply inherited. The underlying algebras are given by 
the corresponding representations of (|IJ) . 

Evidently, fixing M~ = 0,M + = 1. (Q) leads to the well known quantum algebra U q (su(2)) Now the simplest 
representation S = a derives the integrable XX Z spin chain ]l3| , while ^ with the corresponding reduction of g{u) 
p0| yields the lattice sine- Gordon model j2^] with its Lax operator obtained from (|^) with all c's = f . 

An asymmetric choice of central elements: c\ 2 = f, c\ = -4r = —iq,q — e la along with the mapping S + = 

cA, S~ = cA\ S 3 = —N, c = (cot a) 2 brings (|f|) directly to the well known (/-oscillator algebra p3| , p4| and 
simplifies g 2 {u) = [— 2u] q in (^|). Therefore using the inter-bosonic map |25|] one gets a bosonic realization for the 
g-oscillator |2^]. This realization in turn constructs easily from (||) the Lax operator, which coincides exactly with the 
discrete version of the quantum derivative nonlinear Schrddinger equation (QDNLS) Jl^ |. The QDNLS was shown 
to be related to the interacting bose gas with derivative 5-function potential pq |. Fusing two such models one can 
further create an integrable massive Thirring model described in Q . 

Having the freedom of choosing trivial eigenvalues for the central elements: c{ = ci, = with other c's— 1 we 
obtain another deformed Lie algebra [S + , S~] = e isina ■ This can be realized again by (||) with the related expression 
for g{u) [20| , using which the Lax operator is obtained from (Q). The model thus results is no other than the discrete 
quantum Liouville model ^7j. Note that the present case: M ± = ±-\/±l may be achieved even with c{ 7^ giving 
the same algebra and hence the same realization. However, the Lax operator which depends explicitly on c's gets 
changed reducing (Q) to another nontrivial structure. This is an interesting possibility of constructing different useful 
Lax operators for the same model, in a systematic way. For example, the present construction of the second Liouville 
Lax operator recovers that of [Q, invented for its Bethe ansatz solution. 

In a similar way the particular case M = can be achieved with different sets of choices: with all c's = except 
«)c+ = I , ii)cf = ±1, iii)cf = I, all of which lead to the same algebra 

[S+,S-] = 0, [S" 3 ,^] = ±S ± . (4) 

However, they may generate different Lax operators from (||), which might even correspond to different models, though 
with the same underlying algebra. In particular, case i) leads to the light-cone SG model, while ii) and iii) give two 
different Lax operators found in |28| and |2j| for the same relativistic Toda chain. Since here we get g(u) =const. , 
interchanging u — > ~ip,p — > —iu, (||) yields simply S 3 — —ip, S ± = ae TU generating discrete-time or relativistic 
quantum Toda chain. 

Remarkably, all the descendant models listed above have the same trigonometric i?-matrix inherited from the 
ancestor model and similar is true for its rational form, as we will see below. This unveils the mystery why a wide 
range of models found to share the same ii-matrices. The L-operators and the underlying algebras however become 
different, being various reductions of the ancestor Lax operator (pi) and the ancestor algebra (pi). 
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We consider now the undeformed a — > limit of the proposed algebra ([j]) . It is evident that for the limits to be 
finite the central elements must also be a dependent. A consistent procedure leads to — > is ± , M + — > — m + , M~ — > 
—am - , £ — > 1 + iaA giving the algebraic relations 



[s + ,s>2mV+m-, [s J ,s ± ]=±s ± (5) 

with m + = c\c2,ra~ = c\c\ + c5c| as the new central elements. It is again not a Lie but a QdA, since multiplicative 
operators rri^ can not be removed in general due their allowed zero eigenvalues. (^|) exhibits also noncocommutative 
feature Jl^] unusual for an undeformed algebra. (||) at this limit reduces to 



c?(A 



s 



3\ 



MA)= \\ , * , , (6) 



c^(A- S . 

while i?-matrix is converted into its rational form, well known for the NLS model [pj . Therefore the integrable systems 
associated with algebra (||) and generated by ancestor model (||) would belong to the rational class all sharing the 
same rational i?-matrix. 

It is interesting to find that the bosonic representation (^), using the undeformed limit go(u) [pc|l and the inter- 
bosonic map [^5| reduces into a generalized Holstein-Primakov transformation (HPT) 

s 3 = s-N, s+ = g {N)ip, s~=^go(N), gl(N) = m~ + m + (2s — N), N = ip^ip. (7) 

It can be checked to be an exact realization of (||), associated with the Lax operator (^). This would serve therefore 
as an ancestor model of rational class and represent an integrable generalized lattice NLS model. 

For the choice m + — l,m~ ~ 0, (||) leads clearly to the standard su(2) and for spin ^ representation recovers the 
XXX spin chain [0. On the other hand the general form (Q) simplifies to standard HPT and (^|) reproduces the 



lattice NLS model |22 



The complementary choice m + = 0,m~ = 1, reduces (|^) to a non semi-simple algebra and gives go(N) = 1. This 
induces a direct realization through oscillator algebra: s + = ip,s~ — tp^s 3 = s — N and corresponds to another 
simple lattice NLS model |}(|. Remarkably, further trivial choice ra~ = gives again algebra (||) and therefore the 
same realization found before for the relativistic case can be used, but now for the nonrelativistic Toda chain Q . The 
associated Lax operator should however be obtained from (^j) along with the rational i?-matrix. 

It should be noted that a bosonic realization of general Lax operators like (||) and (||) can be found also in some 
earlier works p0| , pl| . Apart from the discrete models obtained above, one can construct a family of quantum field 
models starting from their lattice versions. Scaling first the operators like pj, Uj, cj, ipj, consistently by lattice spacing 
A and taking the continuum limit A — > one gets pj —> p(x),ipj — > ip( x ) e ^ c - The Lax operator C(x, X) for the 
continuum model then obtained from its discrete counterpart as Lj(X) —* I + iAC(x). The associated i?-matrix 
however remains the same since it does not contain A. Thus integrable field models like sine-Gordon, Liouville, NLS 
or the derivative NLS models are obtained from their discrete variants constructed above. 

It is possible to build further a new class of models, that may be considered as the inhomogeneous versions of 
the above integrable models. The idea of such construction is to take locally different representations for the central 
elements, i.e. instead of taking their fixed eigenvalues one should consider them to be site dependent functions. This 
simply means that in the expressions of g(uj) |2(|, should be replaced by and consequently in Lax operator 
(^|) all c's should be changed to c^-s. Thus in lattice models the values of central elements may vary arbitrarily at 
different lattice points j including zeros. This would naturally lead to inhomogeneous lattice models. However since 
the algebra remains the same they answer to the same quantum i?-matrices. Physically such inhomogeneities may be 
interpreted as impurities, varying external fields, incommensuration etc. 

Notice that in the sine-Gordon model unlike its coupling constant the mass parameter enters through the Casimir 
operator of the underlying algebra. Therefore taking Mf = —(Amj) 2 , one can construct a variable mass discrete SG 
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model without spoiling its integrability. In the continuum field limit it would generate a novel sine- Gordon model 
with variable mass m(x) in an external gauge field 6{x). In the simplest case the Hamiltonian of such model would be 
TL = J dx (m(x)(ut) 2 + m~ 1 (x)(u x ) 2 + 8(mo — m(x) cos(2cra))) . Similar models may arise also in physical situations 
II- 

Inhomogeneous lattice NLS model can be obtained by considering site-dependent values for central elements in (^) 
and in the generalized HPT (R), where time dependence can also enter as parameter. As a possible quantum field 
model it would correspond to equations like cylindrical NLS J33| with explicit coordinate dependent coefficients. In 
a similar way inhomogeneous versions of Liouville model, relativistic Toda etc. can be constructed. For example, 
taking c\ — > c°- in nonrelativistic Toda chain we can get a new integrable quantum Toda chain with inhomogeneity 

having the Hamiltonian H = J2j(Pj + ~i) 2 + "o — e Uj ~ Uj+1 . 

Another way of constructing inhomogeneous models is to use different realizations of the general QdA (111) or (H) 
at different lattice sites, depending on the type of i?-matrix. This may lead even to different underlying algebras 
and hence different Lax operators at differing sites opening up possibilities of building various exotic inhomogeneous 
integrable models. Thus in a simple example of impurity XXX spin chain if we replace its standard Lax operator 
at a single impurity site m by a compatible L am = (A + c^Jcr^, the Hamiltonian of the model is modified to H = 

- (T,jjim,m-1 Oj^J+l + ftm-lm+l) > where h m -l m+ l = -(^Ul<i+l+ CT m-l CT m+l)+ CT m-l CT m+l- lt § ives an integrable 
quantum spin chain with a defect, where the coupling constant has changed sign at the impurity site. If attempt is 
made to restore the sign it appears in the boundary condition. 

Thus we have prescribed an unifying scheme for quantum integrable systems, where the models can be generated 
systematically from a single ancestor model with underlying algebra ([!]) . The Lax operators of the descendant models 
are constructed from @ or its q — > 1 limit (^), while the variety of their concrete representations are obtained 
from the same general form (^) at different realizations. The corresponding underlying algebraic structures are 
the allowed reductions of (|l|). The associated quantum i?-matrix however remains the same trigonometric or the 
rational form as inherited from the ancestor model. This fact also reveals a universal character for solving the 
models through algebraic Bcthc ansatz (ABA) [§j5i]|. The characteristic eigenvalue equation for the ABA is given 
by A m (A) = a(X) YYJLi f(^i — ^) + /3(A) YljLi /(A — A 3 -) , where the coefficients a(X) and /3(A) are the only model- 
dependent elements, as being eigenvalues of the pseudovacuum they depend on the concrete form of the Lax operator. 
The main bulk of the expression however is given through functions like /(A) = f^x}) i-e. as the ratio of two elements 
of the i?-matrix and hence is universal for all models belonging to the same class. Therefore all integrable models 
solvable through ABA can be given by almost a universal equation based on a general model. 

The author thanks Prof. Y. Ng for his kind hospitality and encouragement at N. Carolina Univ. where this work 
was completed and Prof. Alexios Polychronakos for fruitful comments. 
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